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Abstract 

We investigate tiie optimal strategy over a finite time liorizon for a portfolio of stock and bond and a 
derivative in an multiplicative Markovian market model with transaction costs (friction). The optimiza- 
tion problem is solved by a Hamilton-Jacobi-Bellman equation, which by the verification theorem has 
well-behaved solutions if certain conditions on a potential are satisfied. In the case at hand, these con- 
ditions simply imply arbitrage-free ("Black-Scholes") pricing of the derivative. While pricing is hence 
not changed by friction allow a portfolio to fluctuate around a delta hedge. In the limit of weak friction, 
we determine the optimal control to essentially be of two parts: a strong control, which tries to bring the 
stock-and-derivative portfolio towards a Black-Scholes delta hedge; and a weak control, which moves 
the portfolio by adding or subtracting a Black-Scholes hedge. For simplicity we assume growth-optimal 
investment criteria and quadratic friction. 

1 Introduction 

An idealised model of investment is a sequence of gambles where an investor at each time step decides 
■ ■ ■ if to re-balance her investments, and, if so, by how much. The game is multiplicative if the pay-off is 
proportional to capital, and Markov if the new capital and new position only depend on the previous state 
and the action taken then. In two previous contributions we computed the strategy an investor should 
use to maximize the growth rate of her wealth, in the presence of transaction costs, if she can invest in stock 
and bonds. In this paper we extend the investment possibilities to also include a derivative security, e.g. an 
option on the stock. 

Asset allocation optimization in the presence of transaction costs has a long and distinguished history in 
finance. The main mathematical tool is the Hamilton-Jacobi-Bellman equation, introduced in the friction- 
less case by Merton [.22. .23J . and with friction by Constantinides [10|. A pedagogical introduction to the 



1 



application of the Hamilton-Jacobi-Bellman equation to financial problems can be found in Our first 
paper, [3 1, overlaps with the subsequent work of Atkinson, Pliska and Wilmott [1 J, where also a perturbative 
expansion around the friction-less limit is performed. In [|41 we treated also optimization over a finite-time 
horizon. 

Turning to option pricing and hedging in the presence of transaction costs, in several contributions, the 
Hamilton-Jacobi-Bellman approach was not used. A classical paper in this line of research, leading to a 
modified option price, is that of Leland |21 1, another one, also leading to a modified option price, is that of 
Bouchaud and Somette [9|. Boyle and Vorst ^ looked at a portfolio replicating an option with transaction 
costs. This is a different setting to that of optimal control and an Hamilton-Jacobi-Bellman approach. 
Avellaneda and Paras [2J did consider Hamilton-Jacobi-Bellman , but in the limit of large transaction costs. 
Other notable contributions are those of Edhirisinghe, Naik and Uppal iQ, Davis, Panas and Zariphopolou 
ll^j, and Bensaid, Lesne, Pages and Scheinkman [5|. 

A recent contribution considering a portfolio including stock and derivatives in the Hamilton-Jacobi- 
Bellman approach is the paper of Constantinides and Zariphopolou 1 11 1. However, these authors consider 
derivatives which can be traded only once (section 3 of that paper), while we address the situation where 
derivatives are traded continuously. Constantinides and Zariphopolou show that expected utility is increased 
by including derivatives, and that derivative prices must obey certain bounds in their model. We note that 
Constantinides and Zariphopolou 1 1 1 1 remark that the bounds they derive would be tighter if derivatives 
would be traded more than once. 

In our treatment, the Hamilton-Jacobi-Bellman equation for the stock, bond and derivative problem 
leads to time-dependent controls, as in the simpler case of |4|. A feature which appears when stock and 
derivative are both traded continuously is that the optimization problem in general ill-defined, unless the 
price process obeys a solvability criterion, known as the verification theorem in the mathematical theory of 
controlled stochastic processes [ 1 6 1 . In Hamilton-Jacobi-Bellman language, the value function is potentially 
unbounded, because the number of variables that can be controlled (the positions in stock and derivative) is 
larger than the number of independent noise sources. In the case at hand, the potential is the expected utility 
of the portfolio as a function of the fraction of wealth invested in stock and derivative, and the conditions 
simply imply arbitrage-free ("Black-Scholes") pricing of the derivative. 

While we hence find that pricing does not depend on market friction, the optimal investment strategy 
does. Qualitatively speaking, we determine the optimal control to be of two parts: a strong control, which 
tries to bring the stock-and-derivative portfolio close to a Black-Scholes delta hedge; and a weak control, 
which moves the portfolio by adding or subtracting a Black-Scholes hedge. The rationale for the presence 
of the weak control is that the strong control acts to oppose the underlying diffusion of the portfolio, in the 
direction normal to the delta hedge. The larger that diffusion, the higher will be the average friction costs, 
per unit time. It is therefore advantageous to invest as much in the delta hedge to make the diffusion in the 
normal direction as small as possible. As in fm we find that introducing derivatives increases expected 
utility. Essentially, in the small friction limit, expected transaction costs decrease, and dependence upon 
friction parameter is pushed to higher order. 

A technical contribution in this paper, analogous to ID SI (without derivatives), is that we introduce a 
multi-scale expansion around the friction-less limit. Since, however, we have two independent variables 
under control (the stock and the derivative), we can have different scales in two different directions. In 
fact, we will show that in the weak-noise limit there is a fast control direction, and a slow control direction. 
The fast control strives to bring the portfolio to an optimal stock portfolio plus a Black-Scholes hedge. 
Financially, this means that an optimal investment strategy is to hold some amount in stock, and then some 
number of fully hedged derivatives. That number is however controlled on a longer time scale, by the 
slow control. Two limit cases are of interest. First, far from expiry a derivative is not much different from 
stock, and the situation is similar to only investing in stock and bonds. Second we can also deal with the 
situation close to expiry. There, the best strategy turns out to be to hold little funds proportional to the 
Black-Scholes hedge, i.e. to avoid derivatives in the optimal strategy. This concurs with the practice of 
closing out positions in derivatives before expiry. 

For simplicity we work in this paper with quadratic friction. These can be motivated as an effective 
description of market impact (see e.g. |15|). The reader is referred to |4| for details. Linear friction costs, 
arguably more realistic, lead to free boundary problems in the Hamilton-Jacobi-Bellman formalism, which 
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are considerably harder, from the analytical and numerical point of view. 

For simplicity we will also furthermore assume throughout that an investor strives to optimize expected 
growth of capital, which in a multiplicative market model means logarithmic utility. Growth optimal strate- 
gies were first introduced by Kelly in the context of information theory |20|. Growth-optimal strategies 
have the well-known property of eventually, for long times, outperforming any other strategy with proba- 
bility one ([ 18| and references therein), but do not maximize vanilla- flavored utility functions, see e.g. \ 131 . 
In the present context, logarithmic utilities should merely be looked upon a definite and convenient choice, 
which brings some mathematical tidiness. 

The paper is organised as follows. In section|2lwe state the model (without derivative), and the controls 
we consider. We state the optimisation problem in the framework of the Hamilton-Jacobi-Bellman equation. 
In sectionlUwe show that the non-linear Hamilton-Jacobi-Bellman equation governing the dynamics, in our 
example, is solvable in the small transaction costs limit by means of a multi-scale perturbation theory 
(see for example |8|, or |17|, chapter 9). This is the main technical result of the paper, and reduces the 
non-linearity to a normal form. All higher order corrections can be computed from ancillary linear non- 
homogeneous equations. In section |6l we solve analytically the normal form of the non-linearity. The 
approximation turns out to be very accurate for realistic values of the parameters in the model. The last 
section is devoted to a discussion of the results. 



2 Bond, stock and derivatives 

In this section we define notation, and state the problem. The market consists of a risk-less security (bond, 
or bank account) and a risky security (stock). By a change of numeraire we take the price of the risk-free 
security to be constant in time. The stock price is taken the standard log-normal process: 

i/jt+dt = i't{l + pdt + adBt) (1) 

Here dBt denotes the stochastic differential and fi and a are positive constants. Nothing in the following 
analysis would essentially change if /i and a would be functions of t and St, as long as the market is still 
complete, see e.g. 1 19|. Consider now first a portfolio in only stock and bond. The control variable then is 
the fraction of wealth an investor has invested in stock: 

TT/Stocks 

This variable changes both in result to market fluctuation, i.e. ([T]), and re-hedging. We assume that a control 
can be executed of the form 

= fdt (3) 

and doing so carries a cost 

^p^tradingcost. ^ 

where F is a semi-positive definite functional of the stochastic control. The form of F models the transac- 
tion costs. 

The coupled stochastic differential equations of W and p are then 

dWt = W[{mdt + aptdBt-jF{f)dt)] (5) 
dpt = [f + Pt{l - Pt){p - (^^Pt)] dt + [apt{l - Pt)] dBt + lPtF{f)dt (6) 

For a derivation of these equations, see ||U. The time-dependent growth-optimization problem, of a stock 
and bond portfolio, from time t to some final time T, is simply to choose the control / such that the 
expected value of log ^ is maximized. By a change of variable, this is equivalent to maximizing the 
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expectation value, over the controlled diffusion process, of a potential (utility function) depending on p and 
/. Implicitly, we assume unbounded borrowing at the risk-less rate, and no restrictions on going short. 
These are not in fact serious limitations, because the optimal solution, with transaction costs, is to hold the 
fraction invested in stock close the optimal value of ^, see |3|, which is finite. In the main body of this 
paper, we will use quadratic friction costs, i.e. F(f) = p, for a discussion of linear friction costs, see [l3||. 

Let us define a derivative security as a third investment possibility, the price of which, C{'ipt, i), only 
depends on the moment of time t and the price of stock. The price dynamics of the derivative is 



dC 

~c 



1 
c 



dtC + fiiptd,ptC + 



C 



jjddt + ad dBt 



(7) 



where we for later convenience introduce amplitudes jid and ad- Both are of course functions of t and ipf 
Let now as before the fraction invested in stock be p with control /, and the fraction invested in derivative 
Tj, with control fd- Exercising either of the controls in a time interval dt carries a cost fd)Wtdt. 
The coupled equations for wealth, p and rj are then 

dW = W[{fip + fidV)dt+{ap + adr])dBt-J'{f,fd)dt] (8) 

dp = [f + a + pT{fJd)]dt + hdBt (9) 

dri = [fd + ad + V -^(/, fd)] dt + bd dBt (10) 
where the functions in the drift terms are 

a:= fip- p [pp + fidV - {crp + cfdV - o") (crp + adif])] (11) 

ad = PdV-V [yUP + fJ-dV - {<^P + o-dV - cfd) {crp + adr])] (12) 
and the functions in the diffusive terms are 

b = ap- p {ap + adT]) (13) 

hd = (^drj - T] {ap + adT]) (14) 

With analytic transaction costs we have 

HfJd)=l\f\' + ld\fd\' (15) 

with two friction parameters 7 and 7^. We now state the problem we want to solve. The expected logarith- 
mic growth rate is 



log 



T 



Wt 



X{x,y,p,t;T) := E 

pt=x;rit=y;ipt=p 

In consequence the logarithmic growth is the expected value of the utility function 

{ap + adVy 



U = p,p + pdV 



(16) 

(17) 
to reach 

it was at (x, y, p) at time t: 

\{x,y,p,t;T) = dt' j U{z)P{x',y',p',t'\x,y,p,t)dx' dy' dp' (18) 

Note in view of (HJ) the probability density factorizes to 

Pix',y',p',t'\x,y,p,t) = Ppn{x',y',t'\x,y,p,t)P^{p',t'\p,t) 

Furthermore, the probability density is in general non-autonomous as pd and ad may depend explicitly upon 
the time variable. The problem is now to find controls / and fd that maximize the logarithmic growth. 



over the probability density P{x' ,y' ,p' ,t'\x,y,p,t) is the probability of the process (pt,r]t,tpt) 
point (x', y', p') at time t', given it was at (x, y, p) at time t: 
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3 The verification principle and Black-Scholes 



It is useful to first discuss the friction-less case. We will then just reproduce standard elementary results 
in finance, but in a formulation convenient for the following discussion. Without transaction costs, the 
speculator is free to rehedge continuously. In such a case the optimisation problem is equivalent to finding 
the supremum, at any instance of time, of the instantaneous growth rate 



V = IJ.P + fXdV - 



(19) 



Equation ([T9l) is a degenerate quadratic functional of the fraction in stocks and derivatives. The Hessian of 



H 



-a 



has a zero eigenvalue along the marginal subspace 

a p + adf] = 

The second eigenvalue is negative, hg = — (cr^ + a^), associated to the stable subspace 

CTdP - cr?7 = 

We now make a change of variables 



peio + r] eoi 



(20) 



(21) 



(22) 



(23) 



where (eio, Coi) is the canonical basis of and e^) is an orthonormal basis formed by the unit vectors 
respectively spanning the marginal and stable subspaces of the Hessian matrix H: 



—a 



VWs 



a 



(24) 



The variable ( along the stable eigenspace describes a portfolio in which the investment in derivatives is 
weighted by the ratio of the volatilities 



C = P+—V 

a 



The utility function reads in these new variables 



V = IJ'C+ [Pd CTd 

a 



(25) 



(26) 



This growth rate is a convex function if and only if the second term vanishes. This can happen if either of 
its two factors are zero. The first possibility gives the following solvability condition: 



a 



fi pd = 

o-d 



df.C 



C - 



(27) 



holding for every t' G [t, T] and in particular for f equal to t: 

dtC+^d'^C = Q (28) 

This is the of course simply Black-Scholes equation at zero interest rate. The second possibility is that the 
linear combination 

adC-a^ = (29) 

vanishes, which simply means that the fraction invested in derivatives is zero. Optimisation can then be 
carried out along the stable manifold. The utility 

V = fiC-^ (30) 



has a maximum for 



c = 4 (31) 



If nothing is invested in derivative {r] = 0) the fraction invested in stock (p = ^) is the same as the optimal 
investment fraction in the stock-and-bond problem. When dTTl holds true any dynamics along the marginal 
subspace does not produce any gain or loss. That means we can invest p W (in value) in stock and —-^pW 

(in value) in derivative, for any p. Expressed in stock price tpt and numbers of stock, n^, the value invested 
in stock is n^ipt- The value invested in derivative is hence —^n^tpt, but also n^C, if is the numbers of 
derivative. The number of stock per derivative is hence 

— = = -a^C (32) 

Ud a tj) 

The portfolio along the marginal subspace is hence a simply Black-Scholes delta hedge 

/\:=d^C (33) 

following the standard financial notation. 



4 Hamilton-Bellman- Jacobi problem for analytic transaction costs 

The use of analytic transaction costs renders the Hamilton- Jacobi-Bellman problem simpler to study. In the 
frame of reference fixed by the eigenvectors of the Hessian of the utility function, the stochastic dynamics 
is governed by the system of stochastic differential equations 



dW = W [pC dt + aC dBt] - WT (/^ /^) dt (34) 

c?C= [/c + ac + C^(/c,^)]^^ + &C^^* (35) 

d-d = [U + a^ + '&J' (/c U)] dt + h dBt (36) 

dad = (y H dt + a K dBt (37) 

dijj = pi/j dt + a ijj dBt (38) 
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The drift and diffusion fields in these coordinates are 



C(l - C) + crlCi'Jd - 'JQ + - <Jd) 



adC - a'd 



[H + K{(Ja-(TC)] 



aH{l-C) + <ylC{aa-aQ + aaad{a-aa) , ^dC"^ ^ 

r (J — — A 



adder - ad) 



(39) 
(40) 



a' + a, 



2 +^[<^d ^2t?2 - C 



+a 
b^^a 



[H + ail-C)K] 



+ a — ^ —K 



a' + a'. 



(41) 
(42) 



while the time change of ad is expressed in terms of two new amplitudes: 



H 



dtad + d^ad H T^d^ad 



K^^d^ad (43) 
By using the Black and Scholes equation, H can be expressed in terms of K and the other parameters: 

= + ^) (44) 

The dynamic programming equation for quadratic friction 

Hk , U) = 77^ (/c + ^^)' + Id^ (^/c - U)' (45) 



IS 



^^«2x , 



(46) 



The optimal investment strategy is sought by requiring the capital growth A reach as a functional of the 
controls (/^ , /^) a stationary point: 



dT 

+ z—d,X + ( 1 



5/c a/^ 



5/^ 



y a,A = 



(47) 



The stationary point equations admit a unique solution for the stochastic controls 

i.l'^d + Id (^^) dzX + aad{ld - l)dy\ 



2 7 7d (1 _ zd:,X - ydyX) 



^ (7^^ + Id (yf) dyX + aadild - l)d,X ^^^^ 
2 77^^2 (1 - zd,X - ydyX) 

which inserted in the dynamic programming equation yield the the Hamilton-Jacobi-Bellman equation for 
the model 

dtX + ai^d^X + aadyX + ^ipdpX + — ^ — "—^ — — ^ -— — "—^ 

4 a^-f-fd (1 - zd^X - ydyX) 

+^dlX + -^d^X + ^ajA + h^hd^yX + 6cf^p5.pA + hapdy,X 

The dynamics is fully specified by associating to ( B^ the equation for the derivative volatility 

dtad + ^-d,al + ^^Ja, = (50) 

Equation (l49l) contains in principle the solution to our problem. The expected capital growth rate A deter- 
mines the stochastic controls through (l48b . Qualitative and quantitative analyses of the control strategies 
are pursued in sections 15] and l^below. 

The Hamilton-Jacobi-Bellman equation should be supplemented by boundary conditions on A. Ar- 
guably, the most natural would be to additionally assume that the process is confined to some domain, and 
hence that the normal gradient of A vanishes on the boundaries of that domain. For the rest of this paper, 
we will be concerned with a description close to the optimum. We will hence assume that the boundaries 
are far away, that the probability is there small, and that we will not need to further specify the boundary 
conditions. 



5 Qualitative analysis of the Hamilton-Jacobi-Bellman equation 

Qualitatively, one can imagine the following scenario. As the utility changes with the position along the 
stable direction, one may postulate difast control along that direction, striving to bring the position close 
to the marginal. Depending on where one is along the marginal, there would then be different (expected) 
transaction costs in keeping the position in the stable direction close to zero. Hence, all positions along 
the marginal are not equivalent, because they lead to different (expected) transaction costs in the stable 
direction. In fact, we can then posit that the preferred position along the marginal is such that 6^ of (l40b . 
evaluated at 

c = 4 (51) 

is as small as possible. Since this function is linear in the marginal coordinate (see below), we then have 
a prescription for the optimum allocation into a Black-Scholes hedge, as function of time and price of the 
underlying 

bl = D + (cTrfC - (y^ ) YT—2 (52) 
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with 



D:=^^fl-4) (53) 



the diffusion amplitude in the absence of derivative trading. From the definition of K (see appendixlBj) the 
right hand side of (l52b can also be rewritten as 

(54) 

Note that the time variation of the derivative price is usually denoted in the financial literature by the Greek 
letter 0. The relations satisfied by the derivative-© with other commonly used financial indicators as the 
derivative A and T are recalled in appendix |b1 The condition 

bl = (55) 

is enforced by setting 

,^^Eic-(A±^ ,56) 

a 2adtC 

This equation is in one sense the main result of the present work. It is therefore useful to rewrite it directly 
in the original variables, i.e. the fractions invested in stock and derivative: 

o" / r.^^ cTdDC 11 AD 

p = -^—^(^C + ^d^n = 4 - -iirr = 4 + — F 

cr + cr^ ^ 2(^1^ cr a pi. 

^ / DaC AD 

cr^ + a^ 2dtC adpT 

A consequence of these equations is that if cr^ diverges, p* stays finite, while r]* tends to zero (9* diverges 
as cTfi). This happens for European Call Options close to expiry and at-the-money (see appendix, is then 
proportional to the "Greek" A, and inversely proportional to the option price C). Following Black-Scholes 
hedging directly can then lead to large transaction costs, because of "portfolio flapping" (switching between 
the fully hedged and totally unhedged positions in response to small changes in the underlying). We see that 
from the perspective of optimal investment strategies, this problem does not appear, since such portfolios 
do not contain any at-the-money options close to expiry. 

In the opposite limit of a large investment horizon, the derivative volatility tends to the volatility of the 
underlying. Furthermore the inequality 



dtC 



C 



< 1 (58) 



holds true requiring increasingly large investments in derivatives in order to enforce (l55t . In such limit 
the optimal size of the investment in stocks is also seen to diverge. The reason is that the drift (0^ and 
diffusion amplitude ( H^ along the stable direction become independent of the marginal direction if the 
infinite horizon limit is taken for any arbitrary finite value of In particular ( l39b . (l40b become in the 
portfolio variable Q identical to the drift and diffusion amplitude felt in the stock and bond model studied in 
H|. These phenomena indicate a break-down of the argument used to derive (1551) . Qualitatively one expects 
in this second limit the optimal investment strategy to treat the stock and the derivative in a similar manner 
analogous to the one depicted in 
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6 Systematic analysis of the Hamilton- Jacobi-Bellman equation 

In this section we will use a systematic multi-scale analysis to the Hamilton- Jacobi-Bellman equation. We 
will show that the previous qualitative analysis is well founded. Furthermore, we are also able to treat the 
case when the (putative) optimal position 9* diverges, and, more generally, we can compute the control to 
be exerted on C, and 6. Nevertheless, the main interest here is conceptual, that the previous analysis can be 
systematically justified. 

To start with it is convenient to write the logarithmic growth in the form 

A=^(r-t) + (^ (59) 

The first term corresponds to growth in the absence of transaction costs. The intensity of transaction costs 
is then measured by the two adimensional parameters 

e — a (60) 

7 + 7d 

and 

Si = (61) 

7 + 7<i 

In order to construct an asymptotic expansion around the ideal case of no transaction costs it is convenient 
to shift the origin of the coordinates {z, y) to 

Z^C + Z, y^^* + y (62) 

The Hamilton- Jacobi-Bellman equation can be rewritten as 



dt(p + ac^dz^p + a^dyif jipd^^p 



2 



8£ 1 - (C + z)d,^ - {{^* + y)dyip 



8e 1 - (C + z)d,ip - ydy^ ^ ^ 2 ^^^ 2 

z^ 

^-\h^dzykp h(^apdzpip + b^apdyp(fi — = (63) 

By dimensional analysis, one can motivate the following choice of scales in time, stable subspace (z) and 
marginal subspace (y): 

E°°n /'tttt z z y 

^^^(n) ( ~T/2' TsTs' 7174' TTTs' ■ ■ ■ ' 7174' TTTs' ■ ■ ■ ' 7178' • • • '^'-P 

n=0 ^ 

Introducing fast and slow variables 



tl/2 
Zl/4 



t J. t 

' 3/8 • ^378 5 

JTJI ) 2^1/8 '■= JT/I t 

y 

e 



1/8 ) 
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derivatives are rewritten as 



3/4 



1/2 



+ . . . 



+ e 



7/8 



5.,, g 0(0) +5.1/4 '/'(I) 



+ 



(64) 



We introduce "slower" variables in the perturbative expansion for the following reason. A general feature 
of the expansion we are considering is that lower order approximate solutions enter the linear partial dif- 
ferential equations governing higher order ones in the form of non-homogeneous terms. By Fredholm's 
alternative it follows that the perturbative expansion is consistent if and only if these non-homogeneous 
term have no overlap with kernel of the linear differential operator associated to the homogeneous part of 
the equations. According to the standard multiscale method (see for example [8|) slower variables can be 
used to enforce the consistency conditions. The hierarchy of perturbative equations starts with 



r-l/2 . 



1/2 



+ 



;i + 5^)a^ + (1 - 5^)al 



d 



21/4 



bf , 

2 1/4 



0(0) 



" ^1/4 



(65) 



,5/8 . 



,3/4 



L{tij2, 2:1/4, 0o)0(l) + ^5^2^^(;/)(i) = -L{t3/8, Zi/s, 0o)0(O) 



2 y-^i/i-rm J v^fi/8'^(o) j " ^i/4?/i/8 



where bt, b^ are evaluated at z equal and depend parametrically upon i!)* whilst 



(66) 



L{t,z 



dt 2 



dzi 



/4 



d_ 

dz 



(67) 



7 Leading order asymptotics 

Formally the leading order of the perturbative hierarchy of equations coincide with the one of the stock and 
bond market model studied in [3j. Setting 



A 



(68) 
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for z sufficiently small the asymptotic expression of the logarithmic growth of the investor capital is H 



\{z,t- T) 



,,2 h*^\ /^2\ 1/2 2 



2bfe 



oo —2n 



n=0 




(69) 



with denoting the Hermite polynomial of order 2n. The argument of the logarithm can be further 
resummed using the Fourier representation of the generating function of the Hermite polynomials 



HJax) 



9 9 d^fy 



y=ax 



The result is 

A(z,t; T) 



2a2 



— e 



1/2^ 



{T-t)~e 



-pl/2 



2x1/2 ^2 



(70) 



tanh ■ 



'a'^AT-t 
~2 ^ 



6f £ 



In 



1 + e 



-2i 



t2 a T-t 



(71) 



At variance with |T1 the diffusion coefficient h*^ in dTTT) depends for any finite investment horizon upon 6*. 

The logarithmic growth A attains a maximum for z equal to zero corresponding to the optimal portfolio in 
the absence of transaction costs. The value of this maximum defines the growth rate of the investor capital. 
It is straightforward to verify that the conditions (1551) . ( 15^ specify the supremum for the growth rate of the 
investor capital. The overall logarithmic growth becomes in such a limit 



\{z,t- T) 



6*^0 



{T-t) 



,1/2 



2\ 1/2 2 

— — tanli 
A 2 



>a^AT-t 
2 ^ 



+ 0(.i+VV(i)) 



(72) 



The qualitative conclusion that can be inferred from (1721 is that the inclusion in the optimal portfolio of a 
derivative product quells the effect of transaction costs from the capital growth rate. 

The mathematical conditions for the validity of the asymptotic expression (l72t of the logarithmic growth 
are determined by (l57t . The corresponding portfolio is well defined close to maturity and for values of the 
underlying price close to the strike price, when the volatility of the derivative price becomes very large. It 
is also worth stressing that the asymptotics ^1% holds true for values of z sufficiently small that the effect 
of the boundary conditions can be neglected: 



2 < 1 



(73) 



The reasoning allowing to derive the asymptotics (l72b from dTTT) breaks down in the large investment horizon 
limit as discussed at the end of section|5l Namely in such a limit the terms proportional to in (l65t vanish, 
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leaving with an equation in the portfolio variable z of the same form of the one describing the investment 
strategy in the absence of derivatives the solution whereof was studied in [4J. 

The analysis of the intermediate dynamical regime between maturity and large horizon requires to take 
into account the boundary conditions associated to the Hamilton- Jacobi-Bellman equation (l49t and is be- 
yond the scope of this paper. 

7.1 Corrections to the leading order 

Inspection of (l66t shows that it is consistent to set 

0(1) = (74) 

with 0(0) independent of the first set of slower variables. Hence the first non-trivial correction to (l72b turns 
out to be of the order 0(5^+^/^0(2) ) as in the case of a market model without derivative products H. 

8 Conclusions 

We have shown that optimum investment strategies in a portfolio of stocks, bond and a derivative can be 
determined by Hamilton-Jacobi-Bellman techniques. Black-Scholes equation appears as a solvability con- 
dition for the problem to be well-founded. Optimal strategies can be described as "fuzzy Black-Scholes": 
if transaction costs are small, optimal portfolios are not far from Black-Scholes delta hedges. 

We believe it of conceptual interest that Black-Scholes pricing emerges as a solvability condition for an 
ensemble of possible investment strategies. Hence, Black-Scholes has been motivated in a weaker setting, 
where there is no replicating portfolio. Second, we have shown that expected transaction costs can be low- 
ered by choosing between investments in both stocks and derivatives, and not only in stock. This is not 
surprising, but the point has not previously been made previously by systematic analysis, to our best knowl- 
edge. We note that the qualitative analysis can be extended to the case of several derivatives on the same 
stock. Although there is a "law of diminishing returns", expected transaction costs can then be lowered 
further. 

Finally we have made explicit the optimal fraction invested in derivatives in terms of the standard finan- 
cial "Greeks". 

A European call option 

The boundary condition associated to Black and Scholes's equation for an European call option is 

C(p, 0) = max {p - p, 0} (A.l) 
with p the exercise (strike) price. The solution at zero discount rate is 

C(p, T-t)= piV(0i) - piV(02) (A.2) 

where 



N{x) = dy 




1 + Erf 



X 

7! 



(A.3) 
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and 



ln(| +^ 



a 



Vr^t (A.4) 



Observing that 



dp dp 



(A.5) 



it is found that for an European call option 



_ apN{<p{) 

pN{<P,)-pN{<P,) ^""-^^ 

Thus, in the large investment horizon limit T — t ] oo and in the limit of underlying prices much larger than 
the strike at maturity date p ^ p the volatility of the derivative product tends to the volatility of felt by the 
underlying. 



B Relation with the "Greeks" 

The sensitivity of the derivative price to variation of the underlying are measured by the "Greeks": a set of 
factor sensitivities used extensively by traders to quantify the exposures of portfolios that contain options. 
In the present case the Greeks of relevance are 



dC ^ _ d^C ^ _ dC 
dp ' dp^ ' dt 



A-—, r:=^, 6:=— (B.l) 



In such a case 



(Td = (B.2) 



C 

From the definition (l43t of the field K and the Black and Scholes equation (l28b it follows 

pA p2^2 



(B.3) 



a J aC 



For an European call options 



ippe '^^(T-t) 



(ln£)2 2 

e1 i^(T-t) 



e = -W— ; (B.4) 



the fields K tends to a distribution when t tends to the maturity date. 
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C Asymptotics of the probability distribution of the stock investment 



The leading order asymptotics to the Hamilton Bellman Jacobi equation can be written as 



dtX + T- i^^xf + T^^lx + 



It describes the evolution equation of the potential of a drift field 



V = OzY 

2e ^ 



advecting the Fokker-Planck equation 



dtP + d,{v p) - -d^ p = {] 



(C.l) 



(C.2) 



(C.3) 



describing within approximation the probability density of the investment in stock. The general solution 
can be written in path integral form 



P{z',T\z,t) = J V[C]6{(:t-z)5{CT-z') e 



A 



with 



A=— ds 
Jt 



+ 27^^^ 



(C.4) 



(C.5) 



The path integral can be performed exactly (see ll24l for details) 

_ (z^+z'^) coBh{a)(T-t)}-2z z' 
Afj?(T-t) Ax(z,t) e ^ sinh{u;{T-t)} 

P{z\T\z,t) = e^^^ ~ — ^^^=^^^^ 

a/2 TT sinh u(T — t) 

having used the boundary condition 



and the notation 



~2e 



(C.6) 



(C.7) 



(C.8) 



The explicit form of x is obtained by irnposing probability conservation over z' . If this latter variable takes 
values on the entire real axis, the result dTTT) given in the main text as leading asymptotic to the full solution 
is recovered. The corresponding form of the probability distribution is 



P{z',T\z,t) 



UJ 



2nu tanhtj(T — t) 
A direct calculation allows to verify that (IC.9I ) satisfies the equality 



(z'-z/coBh{a.(T-t)})2 
g 2 1/ tanhaj(T-t) 



T 



ds / dy 



At 

2a^ 



z" 



P{y,s\z,t) 



(C.9) 



(C.IO) 



with xi^i i) also given by dTTl) as required by the stochastic dynamics underlying the Hamilton- Jacobi- 
Bellman equation. 
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